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1.  Proton-Boron  (p-BII)  Colliding  Beam  Fusion  Reactor 

A  detailed  description  of  our  results  has  been  written  [1]  and  will  be  submitted  for  publication.  (A 
copy  is  attached  to  this  report.) 

A  brief  summary  of  our  results  is  as  follows.  We  find  that,  in  the  proton-boron  colliding  beam 
fusion  reactor,  the  power  which  must  be  supplied  to  maintain  an  optimal  colliding  beam 
configuration  is  estimated  to  be  at  least  5.1  times  greater  than  the  fusion  power.  This  implies  that 
effective  power  conversion  efficiencies  to  electrical  power  in  excess  of  84%  will  be  required. 
Furthermore,  if  the  transverse  collisional  spread  of  the  proton  beam  is  to  be  limited  by  electron 
drag,  the  boron  density  is  constrained  to  have  magnitudes  well  below  the  optimal  value  at  which 
the  fusion  power  is  maximized. 

2.  Low-Aspect-Ratio  Belt  Pinch  Research 

The  belt  pinch  is  being  examined  at  the  IFS  as  a  compact  low-neutron  reactor  that  would  be 
potentially  suitable  for  a  naval  power  source.  The  belt  pinch  is  a  tokamak  that  is  extremely 
vertically  elongated,  well  beyond  the  typical  elongations  of  present-day  tokamaks.  In  particular,  it 
has  the  potential  for  achieving  a  much  higher  plasma  beta  value,  which  would  bring  the  power 
density  into  a  range  where  a  1 00MW  reactor  appears  to  be  quite  possible  for  a  large  naval  vessel. 

There  are  four  main  physics  issues  with  the  belt  pinch  approach  that  need  to  be  researched: 

(a)  MHD  vertical  stability  of  the  plasma 

(b)  MHD  kink  stability  of  the  plasma 

(c)  Confinement 

(d)  Neoclassical  tearing  modes  (bootstrap  current  driven  islands) 

We  have  made  significant  progress  in  areas  (a),  (b)  and  (d),  and  some  promising  initial  results 
have  been  obtained  concerning  (c).  Below  we  describe  this  work. 

(a)  MHD  vertical  stability  of  the  plasma 

In  the  area  of  vertical  stability,  major  modifications  were  made  on  an  existing  resistive  wall 
stability  code  at  the  IFS  that  now  enable  it  to  handle  realistic,  highly  elongated  MHD  analyses. 
These  modifications  were  undertaken  because  investigations  with  the  previous  code  revealed 
significant  inadequacies  when  applied  to  extreme  equilibria  such  as  the  highly  elongated  ones 
presently  under  consideration.  These  modifications  have  been  completed,  and  the  code 
performance  is  good.  Using  the  modified  code,  we  have  found  that  stabilization  of  the  vertical 
instability  is  quite  possible  with  an  appropriate  conducting  shell.  There  is  no  obvious 

technological  difficulty  for  such  a  shell  in  a  D-He^  reactor.  (In  a  D-T  reactor,  the  use  of  liquid 
metals  for  the  conducting  shell  might  be  required  for  acceptable  breeding.) 
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(b)  MHD  kink  stability  of  the  plasma 

Ideal  kink  stability  has  been  examined  with  the  use  of  workhorse  codes  obtained  from  the 
Princeton  Plasma  Physics  Laboratory:  viz.,  the  MHD  stability  code  PEST  and  the  equilibrium 
code  JSOLVER.  Our  results  to  date  indicate  that  acceptable  beta  values  can  be  attained  for  naval 
requirements.  Due  to  code  limitations  on  both  JSOLVER  and  PEST,  however,  some  of  the  most 
promising  equilibrium  shapes  cannot  be  examined.  These  investigations  have  proven  to  be 
exceptionally  difficult,  due  to  these  numerical  limitations.  The  extremely  elongated  equilibria 
under  consideration  require  resolution  at  the  limit  of  the  code  capabilities,  which  necessitates 
extensive  checks  of  convergence.  Getting  results  has  consumed  large  amounts  of  computer  time 
and  operator  time.  Nevertheless,  we  have  found  that  ideal  kink  stability  can  be  obtained  with 
plasma  beta  values  of  18  -  24  %  for  elongations  of  4  to  5.  This  would  enable  a  100  MW  reactor 

(including  shielding  for  the  superconducting  magnets)  for  D-He^  (or  D-T)  with  a  plasma  major 
radius  of  about  3  meters  and  a  height  about  twice  that. 

We  believe  that  it  is  possible  to  increase  the  beta  considerably  (probably  up  to  nearly  40%),  but 
different  codes  would  have  to  be  used  to  confirm  this  point.  We  communicated  with  General 
Atomics,  and  they  agreed  to  provide  us  with  the  equilibrium  and  stability  codes  TOQ  and  GATO 
for  this  purpose.  We  have  already  obtained  and  modified  TOQ,  and  we  have  verified  that  it  can 
produce  the  more  attractive  plasma  shapes  that  are  desired.  Next,  stability  analysis  with  the  use  of 
GATO  (or  possibly  the  DCON  code)  will  be  performed  in  order  to  determine  the  degree  of 
improvement  in  plasma  parameters  that  is  possible.  (Some  of  this  work  will  be  done  with  some 
collaboration  with  Princeton  computational  physicists  J.  Manickam  and  J.  Menard.) 

Feedback  stabilization  of  resistive  wall  kink  modes  has  not  yet  been  considered  in  detail,  but 
qualitative  features  of  the  ideal  MHD  studies  indicate  that  feedback  stabilization  is  likely 

practical  (particularly  for  D-He^). 

(c)  Confinement 

Concerning  confinement,  some  nonlinear  gyro-kinetic  simulations  of  electron  temperature- 
gradient-driven  instabilities  have  been  performed  (in  collaboration  with  W.  Dorland  at  the 
University  of  Maryland).  The  results  indicate  that  there  is  a  very  strong  favorable  confinement 
scaling  with  elongation,  even  stronger  than  is  indicated  by  empirical  scaling  laws  for  the 
elongation  regime  1  -  2.  These  scaling  laws  indicate  that  confinement  enhancements  of  1.5  -2 

over  H-mode  would  be  required  for  D-He^  ignition  in  a  100MW  belt  pinch  reactor  (and  no 
enhancement  would  be  required  for  a  100  MW  D-T  reactor).  Present  experiments  find  an 
enhancement  of  up  to  2  under  appropriate  conditions,  when  the  E  _  B  velocity  shear  rate  exceeds 
the  growth  rate.  This  condition  is  very  well  satisfied  for  a  100  MW  reactor.  Thus,  adequate 
confinement  appears  promising. 

(d)  Neoclassical  tearing  modes 

Stability  against  neoclassical  tearing  modes  (i.e.,  bootstrap  current  driven  islands)  has  been 
examined  with  the  more  realistic  equilibria  obtained  from  our  modified  version  of  TOQ.  It 
appears  possible  to  exceed  the  stability  margin  by  a  considerable  factor  for  appropriately  shaped 
equilibria.  This  would  eliminate  one  major  source  of  physics  uncertainty  faced  by  proposed  D-T 
ignition  experiments  (such  as  ITER-FEAT). 

Summary 

In  summary,  initial  results  for  use  of  the  belt  pinch  are  quite  encouraging.  No  insurmountable 
obstacles  have  been  found.  The  belt  pinch  requires  large  extrapolations  over  present  experimental 
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results  for  confinement  and  MHD,  but  these  extrapolations  appear  reasonable  based  on  our 
analysis  to  date  [2], 

3.  Dipole  configuration  research 

Recognizing  that  standard  kinetic-theory  treatments  of  plasma  stability  are  not  valid  at  the  high 
beta  values  (typically  greater  than  unity)  in  a  dipole  configuration,  we  have  developed  an  original 
method  for  analyzing  stability  in  this  regime  [3].  Using  this  method,  we  have  completed  an 
analysis  of  the  stability  of  MHD  ballooning  modes  and  drift  modes  in  a  conventional  magnetic 
dipole  configuration.  (We  analyzed  both  a  laboratory  dipole  experiment  and  the  Earth’s  dipolar 
magnetosphere.)  The  results  indicate  that  a  conventional  dipole  plasma  is  stable  if  the  plasma 
pressure  gradient  scale  length  is  comparable  to  magnetic  field  curvature  scale  length. 

We  are  now  studying  the  possibility  of  creating  stable  transport  barriers — for  example,  by  E_B 
shear  stabilization  and  diamagnetic  stabilization — which  could  potentially  lead  to  the  realization 
of  compact  dipoles.  We  intend  to  pursue  these  issues  and  explore  the  concept  of  compact  dipole 
configurations  in  which  the  plasma  has  stable  confinement  close  to  the  current  ring  that  generates 
the  dipole  magnetic  fields  [4], 
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Abstract 


In  the  proton-boron  colliding  beam  fusion  reactor,  the  power  which  must 
be  supplied  to  maintain  an  optimal  colliding  beam  configuration  is  estimated 
to  be  at  least  5.1  times  greater  than  the  fusion  power.  This  implies  that 
effective  power  conversion  efficiencies  to  electrical  power  in  excess  of  84%  will 
be  required.  Furthermore,  if  the  transverse  collisional  spread  of  the  proton 
beam  is  to  be  limited  by  electron  drag,  the  boron  density  is  constrained  to 
have  magnitudes  well  below  the  optimal  value  at  which  the  fusion  power  is 
maximised. 
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I.  INTRODUCTION 


Rostoker,  Binderbauer,  and  Monkhorst  [1]  (RBM)  have  proposed  a  conceptual  fusion 
reactor  design  which  combines  the  environmentally  attractive  neutron-free  property  of  the 
proton-boron  (p-B)  fusion  reaction  with  the  compact  magnetic  field  geometry  of  the  Field 
Reversed  Configuration  (FRC).  The  reaction  products  are  energetic  charged  alpha  particles 
which  can  be  confined  by  the  magnetic  field.  Since  no  neutrons  are  produced  in  the  primary 
reaction,  radiative  activation  of  the  containing  walls  by  particle  bombardment  is  not  a 
major  concern.  This  design  concept  offers  the  vision  of  a  small  fusion  power  source  with 
the  possibility  for  high-efficiency  direct  energy  conversion  into  electricity  and  with  minimal 
radiation  shielding  requirements. 

In  the  FRC,  beams  of  large  orbit  fuel  ions  propagate  across  the  magnetic  field,  with 
enough  electrons  present  for  the  system  to  be  quasi-neutral.  The  energetic  fuel  ions  follow 
confined  orbits  in  an  annular  field-reversed  layer  with  annular  width  small  compared  to  the 
annular  radius,  and  they  carry  a  significant  fraction  of  the  current  required  for  magnetic 
field  reversal. 

Two  modes  of  operation  are  envisaged: 

(1)  Fuel  ions  move  with  essentially  the  same  cross-field  velocities  and  with  effective 
perpendicular  temperatures  in  “betatron”  motion  of  ~  235  KeV  where  the  fusion  cross- 
section  for  a  thermal  plasma  is  close  to  its  maximum  value.  The  beam  energies  of  protons 
and  borons  are  of  the  order  of  300  KeV  and  3.3  MeV  respectively. 

(2)  Cold  beams  of  fuel  ions  move  at  different  cross-field  velocities  with  relative  energy 
of  collision  at  600  KeV.  At  this  “resonant”  energy,  the  fusion  cross-section  is  a  maximum 
which  is  approximately  three  times  larger  than  the  fusion  cross-section  for  a  thermal  plasma 
of  protons  and  borons.  The  perpendicular  temperature  of  the  beams  must  be  substantially 
less  than  140  KeV  while  the  sum  of  the  beam  energies  of  the  fuel  ions  need  not  be  larger 
than  1  MeV. 

The  proposed  fusion  reactor  concept  presents  formidable  challenges  when  compared  to 
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more  conventional  design  concepts  such  as  the  Tokamak.  The  interaction  energies  of  the 
fuel  ions  are  higher:  ~  235  KeV  in  the  case  of  mode  1,  and  ~  600  KeV  in  the  case  of  mode 
2.  The  energy  gain  per  fusion  reaction  is  smaller:  the  energy  output  per  fusion  reaction 
is  ~  8.7 MeV  (H  +  Bn  =  3  He4  +  8.68  MeV),  and  thus  in  mode  1  the  energy  gain  factor 
is  ~  2.4  if  the  energy  invested  in  the  beam  energy  of  the  fuel  ions  is  ~  3.6  MeV,  while  in 
mode  2  the  energy  gain  factor  can  be  no  more  than  ~  14.5  for  a  resonant  colhsional  energy 
of  interaction  of  ~  600  KeV.  Furthermore,  in  mode  2  operation,  an  optimal  colliding  beam 
configuration  must  be  maintained  in  the  face  of  ion  colhsional  effects:  collisional  heating 
and  expansion  of  the  beam  ions,  and  deviations  from  the  resonant  interaction  energy  due  to 
collisional  drag. 

In  our  review  and  evaluation  of  the  physics  issues  involved  in  (p-B)  fusion  reactors,  we 
will  focus  on  mode  2  operation  since  the  reactor  design  concept  based  on  colliding  beams  of 
protons  and  borons  has  been  the  subject  of  considerable  controversy. 

Lampe  and  Manheimer  [2]  (L  and  M)  have  emphasized  that  the  collisional  rate  upb  for 
(p-B)  momentum  exchange  scattering  is  37  times  faster  than  the  fusion  reaction  rate  vp, 
and  that  the  collisional  (p-B)  drag  rate  although  slower  than  vpb  by  a  factor  of  twice  the 
mass  ratio  mp/m,B  is  still  faster  than  uF  by  a  factor  7.  They  discussed  a  variety  of  colhsional 
processes  which  would  destroy  the  colliding  beam  equilibrium  on  time  scales  much  faster 
than  the  fusion  reaction  time.  They  considered  an  explicit  equation  for  the  time  evolution 
of  the  transverse  spread  of  the  proton  beam,  based  on  the  Fokker-Planck  kinetic  equation, 
and  their  calculations  clearly  showed  the  heating  of  the  proton  beam  due  to  p-B  collisions 
in  the  limit  uvb  >  Vpe,  where  is  the  collisional  rate  for  proton-electron  (p-e)  momentum 
transfer.  They  concluded  that  the  required  colliding  beam  equilibrium  “cannot  be  sustained 
for  long  enough  to  provide  fusion  gain.” 

RBM  responded  that  the  criticisms  of  L  and  M  are  relevant  “if  one  considers  only  colli¬ 
sions  and  omits  the  rest  of  the  physics  of  magnetic  confinement  which  L  and  M  proceed  to 
do  in  the  balance  of  their  report.”  As  stated  by  RBM,  the  problem  with  the  calculation  of  L 
and  M  is  that:  “there  is  no  current,  no  magnetic  field,  and  no  inductance,  and  they  consider 
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only  an  initial  value  problem” ;  “not  a  steady-state  problem” ;  “there  are  no  external  sources 
for  the  protons,  and  the  sink  due  to  burning  is  also  not  considered.”  They  suggested  that  it 
may  be  possible  to  maintain  the  fuel  ion  velocities  by  the  injection  of  pulsed  beams. 

The  comments  of  L  and  M  underline  the  difficult  problem  posed  by  particle  collisions. 
The  fusion  reactivity  decreases  if  the  temperatures  of  the  fuel  ion  beams  acquire  a  thermal 
spread  and/or  if  the  interaction  energy  deviates  from  the  resonant  energy.  Efficient  oper¬ 
ation  therefore  requires  that  optimal  conditions  be  maintained  for  times  long  enough  for  a 
significant  number  of  fusion  reactions  to  occur.  The  technological  issues  involved  in  form¬ 
ing  and  maintaining  an  appropriate  colliding  beam  configuration  are  indeed  very  serious. 
Nevertheless,  the  ’’rest  of  the  physics,”  alluded  to  by  RBM,  does  in  principle  allow  for  the 
possibility  of  quasi-stationary  colliding  beam  configurations  with  lifetimes  longer  than  col- 
lisional  times,  and  it  may  be  somewhat  premature  for  L  and  M  to  declare,  without  further 
investigations,  that  “the  concept  is  fundamentally  flawed.” 

The  ion  beam  velocities  can  be  maintained  constant  by  counterbalancing  the  particle  drag 
forces  with  the  forces  due  to  an  inductive  electric  field  and/or  by  transferring  momentum 
to  the  ion  beams  by  some  appropriate  mechanism  (e.g.  injecting  additional  beams).  The 
control  of  the  transverse  velocity  spread  (“temperature”)  of  the  proton  beam  presents  a 
greater  challenge.  However,  the  rate  of  increase  of  the  proton  ’’temperature”  is  proportional 
to  the  boron  density,  and  at  sufficiently  low  boron  density,  the  proton  ’’temperature”  can 
be  limited  by  electron  drag  when  i 'pb  <  Vpe- 

Given  that  quasi-stationary  colliding  beam  configurations  are  in  principle  not  precluded, 
the  essential  issue  becomes  that  of  determining  whether  or  not  there  are  beam  configurations 
which  could  provide  the  basis  for  a  viable  fusion  reactor. 

RBM  have  analyzed  a  set  of  multifluid  equations  obtained  by  taking  moments  of  the 
Fokker- Planck  equation.  The  particle  distribution  functions  were  taken  to  be  drifted  Maxwell 
distributions  where  the  temperature  and  mean  velocity  may  be  time  dependent.  These  equa¬ 
tions  were  used  to  obtain  design  parameters  for  fusion  reactors  based  on  mode  1  operation. 
These  equations  are  however  inadequate  for  the  investigation  of  reactor  designs  based  on 
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mode  2  operation  since  Maxwellian  distributions  are  not  solutions  of  the  lowest  order  Fokker- 
Planck  kinetic  equations  for  colliding  ion  beam  configurations.  The  merits  of  fusion  reactor 
designs  based  on  mode  2  operation  axe  still  being  evaluated  by  RBM. 

In  this  report,  we  limit  the  scope  of  our  discussion  to  two  topics  which  have  a  direct 
bearing  on  the  the  design  of  colliding  beam  fusion  reactors: 

(1)  We  consider  two  quasi-stationary  colliding  beam  scenarios  subject  to  the  constraint 
that  the  relative  (p-B)  energy  £pb  is  equal  to  the  resonant  energy  of  600  KeV,  one  involving 
an  inductive  electric  field  and  the  other  non-inductive  “external”  forces  acting  on  the  fuel 
ions.  From  momentum  balance  we  estimate  the  minimum  power  which  must  be  supplied  to 
sustain  “steady  state.”  We  find  that  the  power  required  is  significantly  smaller  when  non- 
inductive  “external”  forces  are  used  rather  than  inductive  electric  fields.  Nevertheless,  even 
with  non-inductive  “external”  forces,  the  required  power  is  still  quite  large  when  compared 
to  the  fusion  power;  it  is  estimated  that  the  required  power  must  be  at  least  >  5  times  the 
fusion  power.  In  addition,  it  is  a  major  problem  in  itself  to  sustain  highly  efficient  “external” 
forces  acting  on  the  ion  beams,  and  we  have  not  addressed  this  difficulty. 

(2)  We  explore  the  conditions  required  for  electron  drag  to  be  effective  in  limiting  the 
transverse  velocity  spread  of  the  proton  beam.  Electron  drag  is  effective  if  VpB/vpe  <  1.  At 
“steady  state,”  the  transverse  “temperature”  T±  of  the  proton  beam  is  approximately  given 
by  T±/£pb  ~  VpB/vpe-  Since  T±/£pb  is  required  to  be  small,  we  must  have  vpb / iv  <  1,  an 
inequality  which  can  be  satisfied  if  the  electron  temperature  and/or  the  boron  to  electron 
density  ratio  ns/ne  are  sufficiently  small.  However,  the  fusion  power  is  proportional  to  the 
product  of  the  proton  and  boron  densities  ripUB,  and  with  quasi  neutrality  ne  =  np  +  Zub 
(where  Z  =  5),  the  fusion  power  is  a  maximum  when  nB/ne  =  0.T,  hence  it  would  be 
undesirable  to  have  the  density  ratio  n£/ne  of  borons  to  electrons  much  below  10%. 

This  report  is  organized  into  3  sections.  In  Sec.  II,  we  describe  two  quasi-stationary 
colliding  beam  configurations,  and  we  estimate  the  minimum  power  required  to  sustain 
these  quasi-stationary  states.  In  Sec.  Ill,  we  discuss  the  effect  of  electron  drag  in  limiting 
the  increase  in  the  transverse  velocity  spread  induced  by  proton-boron  Coulomb  collisions. 
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In  Sec.  IV,  we  state  our  conclusions. 


II.  POWER  REQUIRED  TO  MAINTAIN  QUASI-STATIONARY  COLLIDING 

BEAM  CONFIGURATION 

A.  Inductive  Electric  Field 

An  inductive  electric  field  accelerating  electrons  and  ions  in  opposite  directions  can  be 
used  to  balance  the  collisional  drag  forces  which  act  to  decrease  the  relative  velocity  between 
the  particle  species. 

Let  us  consider  a  simplified  model  of  a  colliding  beam  configuration  in  which  beams  of 
electrons,  protons  and  borons  are  in  motion  relative  to  each  other  with  velocities  Ve  =  Vez, 
Vp  =  Vpz  and  V B  =  VBz  respectively.  The  direction  of  beam  propagation  is  taken  to  be 
parallel  to  the  z-axis  ( z  is  the  unit  vector  in  the  2-direction). 

The  beam  configuration  is  assumed  to  be  in  steady  state  on  the  collisional  time  scale. 
The  components  of  the  momentum  equation  in  the  2-direction  with  inclusion  of  an  inductive 
electric  field  E  =  Ez  are 

eE  =  TUpUpeiVp  -  Ve)  +  mBvBp{Vp  -  VB)—  (1) 

Tip 

ZeE  =  mpZ2vpe{VB  ~  Ve)  +  ttib^bp{Vb  ~~  Vp)  (2) 

-eE  =  mpVpeiVe  -  Vp )^  +  mpZ\e(Ve  -  VB )—  (3) 

Tig  Tig 

where  Ze  =  5e  is  the  boron  charge,  ma  and  na  are  the  mass  and  density  of  species  “a” 
(a=e,p,B),  Vpe  and  upB  are  the  proton-electron  (p-e)  and  proton-boron  (p-B)  Coulomb  col¬ 
lision  frequencies  respectively,  and  manauab  —  rribnbVba- 

Note  that  these  equations  are  consistent  with  quasineutrality  ne  —  Up  +  ZnB.  FVom 
Eq.  (1),  Eq.  (3),  and  quasi-neutrality,  we  obtain  the  following  equation  relating  Ve,  Vp,  VB: 

v*-v'=(!t§  )w-v*> 

where  the  parameter  S  is  defined  by: 
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(4) 


s  _  rnB  Tie  VBp 
Tflp  Z  Tip  Vpe 

and  is  proportional  to  the  ratio  of  collision  frequencies  i 'Bp/vpe-  Substituting  in  Eq.  (4)  the 
collision  frequencies  given  by  (see  reference  [2]): 

Vpe  =  2.6  x  10-14ne  (10  KeV/Te)3/2  sec-1 
UpB  =  1.09  x  10_12nB  (100  KeV/Spef 2  sec-1 
uBp  =  9.9  x  10“ 14 np  (100  KeV/SpB)3/2  sec-1 


where  Te  (in  kev)  is  the  electron  temperature  and  SpB  is  defined  to  be 

SPB  =  \mp  (Vp  -  VBf  =  600  KeV 


we  obtain: 


5  =  120.0 


mB 


(: 


H) 

Z  Trip  (  £pB  ) 


3/2 


Expressing  the  relative  streaming  velocities  of  ions  to  electrons  {(V^  —  14),  (Vj,  —  14)}  in 
terms  of  the  relative  streaming  velocity  of  the  ions  (Vp  —  VB),  we  obtain: 


vP-V'=  (|i|)  (Vp  -  vB ) 
Vs  -  ve  =  (ii|)  ( vp  -  VB). 


In  order  to  sustain  the  steady  state,  the  inductive  electric  field  acting  on  ions  and  elec¬ 
trons  must  supply  energy  at  a  rate  sufficient  to  replace  the  streaming  energy  transformed 
into  thermal  energy  by  particle  collisions.  This  power  constitutes  circulating  power  in  a 
reactor  that  is  dissipated  into  thermal  energy. 

The  magnitude  of  the  plasma  current  J  is 

J  =  e  (npVp  +  ZnBVB  -  neVe)  =  ene  {  (Vp  ~  Vb) 

and  the  power  PB  dissipated  by  collisions  is  therefore  given  by: 


Pd 


EJ 

2mBnBuBp 


( Z  +  S f  ne  ,  (1  +  S)2Zne 


[(Z-l)2nBZS  (Z  —  l)2  Slip 


+  1 


'pB- 


(5) 
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The  first  and  second  terms  on  the  right  hand  side  of  Eq.  (5)  are  due  to  electron-ion 
collisions,  while  the  third  term  is  due  to  proton-boron  collisions. 

The  fusion  power  Pj  due  to  (p-B)  collisions  is  : 

Pf  =  npnB  (crv)  £f  (6) 

where  (crv)  is  the  fusion  reactivity  and  £j  =  8.68  MeV  is  the  fusion  energy  released  by  a 
fusion  reaction.  For  polarized  fuel  ions  with  energy  of  collision  equal  to  the  “resonant” 
energy  £pB  =  600  KeV,  the  fusion  reaction  rate  is  (see  reference  [2]): 

(< crv )  =  2.0  x  10-15  cm3  sec-1 . 


The  ratio  of  the  power  dissipated  by  collisions  to  the  fusion  power  is  therefore 


Pd  [  ne  ( Z  +  Sl  rteZjl  +  Sl  \ 

Pf  ‘  \nBZS(Z-l)2  npS(Z-l)2  )' 


(7) 


For  a  given  value  of  nB/ne,  the  power  ratio  Pd/Pj  has  a  minimum  at  S  — 
fi+z(z-i)'nB/w!p/2 1  while  for  a  given  value  of  S  the  minimum  occurs  at  nB/ne  = 

Z-\~S 

Z(Z'/2+S)(Zl/2+l)‘ 

The  power  ratio  Pd/Pj  has  an  absolute  minimum  at  nB/ne  =  l/(2Z)  and  S  =  Z1?2,  and 
its  minimum  value  (Pd/Pj) min  is  : 


4  Z1'2  1 

(£1/2  _  1)2  J- 


For  Z  =  5  (borons),  this  ratio  is  equal  to  ~  35. 

In  Fig.  1  we  plot  Pd/Pj  as  a  function  of  the  electron  temperature  Te  for  nB/ne  = 
0.05,0.1,0.15.  It  has  an  absolute  minimum  of  35  for  nB/nt  =  0.1  at  Te  ~  25  KeV.  The 
required  circulating  power,  much  of  which  is  dissipated  into  thermal  energy,  will  therefore 
be  a  large  multiple  of  the  fusion  power. 

The  power  available  for  conversion  into  electrical  power  is  Pj  +  PB.  If  the  effective 
conversion  efficiency  is  r),  net  electrical  power  is  available  only  if  t)(Pj  +  PB)  is  greater  than 
the  power  PB  required  to  maintain  steady  state,  r}(Pj  +  PB)  >  PB.  For  Pd/Pj  —  35,  the 
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conversion  efficiency  must  satisfy  the  inequality  i)  >  35/36  =  0.97.  Thus,  a  steady  state 
established  by  inductive  electric  fields  is  not  attractive  as  the  basis  for  a  fusion  reactor 
producing  net  electrical  power. 


B.  Non-inductive  External  Forces 


In  an  alternative  “steady  state,”  the  ion  streaming  velocities  can  be  held  constant  by 
“external”  forces,  while  the  electrons  have  an  intermediate  mean  velocity  determined  by 
equality  of  the  opposite  proton  and  boron  drag  forces  on  the  electrons.  It  is  here  assumed 
that  there  is  some  appropriate  and  efficient  means  of  transferring  momentum  to  the  proton 
and  boron  beams  to  make  up  the  losses  due  to  collisional  drag. 

Let  Tp  and  TB  represent  the  equivalent  “external”  forces  which  maintain  the  streaming 
velocities  of  the  ions.  The  modified  momentum  equations  are: 


nB 


Tv  =  mpupe(Vp  -  Ve)  +  mBuBp(Vp  -  VB) — 


Tlr, 


Pb  —  mpZ2vpe(VB  —  Ve)  +  mBuBp{VB  —  Vp ) 


0  =  mp ispe(Ve  -  Vp)—  +  mpZ2upe(Ve  -  VB)—. 

Tle  7le 


(8) 

(9) 

(10) 


The  relative  streaming  velocity  of  the  ions  to  electrons  are  now  related  to  Vp  —  VB  by: 

The  magnitude  of  the  plasma  current  J  is 

J  =  en„  (Vt  -  V,)  +  ZenB  (VB  -  V.)  =  %  ~  Vs)  ■ 

The  dissipated  power,  equal  to  the  rate  at  which  work  is  done  by  the  forces  Tp  and  TB, 
is  therefore  given  by: 


Pd  =  TipVpPp  +  nBVBPB 


TtUbUbUbp 


1  + 


Z  rip 


S(np  +  Z2nB) 
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and  the  ratio  of  dissipated  power  to  fusion  power  is 


1  + 


_ ZUe _ 

S(ne  +  Z(Z  -  l)nB) 


(11) 


In  Fig.  2  we  plot  Pd/Pj  as  a  function  of  the  electron  temperature  Te  for  nB/ne  = 
0.01,0.1,0.15. 

The  dissipated  power  is  lower  than  that  for  a  steady  state  with  inductive  electric  fields. 
For  density  ratio  nB/ne  =  0.1  and  electron  temperature  Te  =  25KeV,  the  dissipated  power 
is  a  factor  of  ~  9  times  the  fusion  power.  The  power  ratio  Pd/P /  decreases  with  increasing 
electron  temperature  but  increases  with  decreasing  density  ratio  nB/ne. 

In  the  limit  S  >  T^z^TjnJ7ne  (or  ^  >  42. 6/(1  +  20nB/ne)2/3  KeV),  the  dissipated 
power  tends  to  its  minimum  value  and  it  is  then  due  essentially  to  proton-boron  collisions. 
Nevertheless,  it  is  at  least  ~  5  times  the  fusion  power  (see  also  reference  [2]),  corresponding 
to  a  required  effective  efficiency  of  conversion  to  electrical  power  of  7?  =  5/6  =  0.83. 

In  Appendix  A  we  establish  that  the  inclusion  of  an  additional  inductive  electric  field 
always  results  in  an  increase  in  the  dissipated  power.  It  is  preferable  therefore  to  maintain 
the  ion  velocities  only  by  non-inductive  “external”  forces  in  order  to  reduce  the  dissipated 
power  to  the  lowest  level  possible. 


III.  LIMITATION  OF  TRANSVERSE  VELOCITY  SPREAD  OF  PROTON  BEAM 

BY  ELECTRON  DRAG 

In  addition  to  maintaining  the  energy  of  collision  of  the  fuel  ion  beams  close  to  the 
resonant  energy,  it  is  essential  to  control  the  increase  in  the  transverse  velocity  spread  of  the 
ion  beams.  Otherwise,  the  annular  width  of  the  field-reversed  layer  will  increase  until  the 
particle  trajectories  intersect  the  nearby  walls  and/or  particles  are  lost  axially  through  the 
x-points  at  the  ends  of  the  FRC. 

In  this  section,  we  investigate  the  effect  of  electron  drag  on  the  transverse  velocity  spread 
of  the  proton  beam  induced  by  (p-B)  Coulomb  collisions.  We  first  consider  time  dependent 
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solutions  of  the  Fokker-Planck  equation  (with  and  without  an  inductive  electric  field)  for 
the  distribution  function  of  a  proton  beam  which  is  injected  into  a  plasma  of  boron  ions 
and  electrons.  We  establish  that  the  velocity  spread  of  the  proton  beam  remains  small  if 
the  (p-B)  collision  frequency  is  less  than  the  (p-e)  collision  frequency,  upb/ vpe  <  1.  We  then 
derive  a  beam  envelope  equation  to  describe  the  time  dependence  of  the  mean  square  width 
of  the  proton  beam  by  taking  appropriate  moments  of  the  Fokker-Planck  equation.  We 
confirm  that  essentially  the  same  condition  is  required  for  small  transverse  velocity  spread 
of  the  proton  beam  when  a  field  reversed  magnetic  field  and  a  transverse  electrostatic  field 
are  included  in  the  analysis. 

We  describe  the  time  evolution  of  the  proton  beam  distribution  function  Fp  in  the  boron 
frame  of  reference  where  the  boron  ions  are  stationary.  The  boron  ions  have  zero  tempera¬ 
ture  while  the  electron  distribution  function  is  approximated  by  a  drifting  Maxwellian  with 
temperature  Te  and  mean  velocity  Ve  =  Vez  which  is  taken  to  be  parallel  to  the  2-axis. 
The  time  evolution  of  Fp(v,  t )  is  determined  by  the  simplified  Fokker-Planck  equation  (see 
Appendix  B): 


f)F  p  f)F 

m+^tpE'l£  =  Cpe{Fp' Fe)  +  Cps{Fp ’  Fb)  +  5(v’ 

where  the  collision  operators  Cpe  and  Cpb  are  approximated  by  : 


(12) 


Cp*  =  -  K.')) 


J pB 


=  VpB  Vq  | 


mTI 


d_p  Va  ,  mB  1  9  1(K  vav0\dFp\ 
(mp  +  mB)  dva  pv3  (mp  +  mB)  2  dva  v  a0  v2  dvp  J  ' 


We  include  a  source  term  S(v,t)  to  model  the  injection  of  the  proton  beam  in  the  z- 
direction  and  an  inductive  electric  field  E  =  Ez  to  balance  the  collisional  drag  forces  on 
the  protons.  We  neglect  proton-proton  collisions.  The  subscripts  a,p  of  vQ,  v$  denote  the 
cartesian  components  x,  y,  z  of  the  velocity  vector  v.  A  stun  over  repeated  indices  a,  (3  is 
implied,  vq  is  the  characteristic  speed  of  the  proton  beam. 

We  introduce  spherical  coordinates  v,  6 ,  <f>  where  vx  =  v  sin  6  cos  (j),  vy  = 
v  sin  6  sin  4>,  vz  =  v  cos  9,  and  we  assume  cylindrical  symmetry  about  the  2-axis.  The 


11 


Fokker-Planck  equation  for  Fp  can  then  be  expressed  in  the  form  : 


dFp  ld3 
Tr  ~  u?  duu  Fp~ 


(  eE 
\TripUpeV  0 


+9  {^l;“2cos^-d^lsin2^} 


VVB 


TUB 


upe  (mp  +  mB) 


{mp  1 

mB  u‘ 


dFn 


+ 


2  du  2u3  sin  6  dd 


dd 


) 


^  sin0^£  \  +S(u,0,t) 


(13) 

(14) 


where  r  =  u  =  v/vq. 

We  consider  the  injected  proton  beam  to  be  monoenergetic  and  sharply  peaked  in  the 
direction  of  propagation,  and  we  model  the  source  term  S  by  : 

S  =  .So©(T)d(it  —  Uj)2aoexp 

where  the  injection  speed  is  given  by  Ui,  the  parameter  a0  1  determines  the  initial  angular 
spread,  and  S  is  finite  only  for  values  of  0  ^  1.  0(r)  is  a  step  function  of  the  time  variable 


r. 


We  are  interested  in  solutions  for  which  the  proton  distribution  function  remains  sharply 
peaked  in  the  direction  of  propagation,  and  we  simplify  the  Fokker-Planck  equation  further 
by  making  the  substitution  sin  6  — *  9,  cos  6  — *  1  —  82/ 2. 


A.  Zero  Inductive  Electric  Field  E  —  0  and  Ve  =  0 

In  the  case  where  there  is  no  inductive  electric  field  ( E  =  0)  and  the  electron  mean 
velocity  Ve  =  0,  we  have: 


—  =  Z27fMu)Fp  +  2^e§e6^  +  S°e(TWu  ~  Ui)2aoexp(-ao02/2 )  (15) 


dr  u2  du ' 


q{u)  =  u3  +  emA 


_  mp 


A  = 


mB 

'pB  1 
^pe  (1  "b  ^m) 


The  solution  of  this  equation  is  [4]  : 


Fp(u,t)  =  Fp(u,t)2cl(u)  exp{-a(u)92 /2} 


(16) 
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where 


Fp(u,  r)  =  3^  {©(u  -  u{ui,r))  -  Q(u  -  it*)} 

u3  +  emA 

(  1 _  a0 

I1  +  ft  loS  {(!  +  ZmA/u*)  /  (1  +  Cmi4/t^)}] 

and  u(ui,t)  is  the  solution  of  the  characteristic  equation  ^  ^  with  initial  conditions 

u  =  U{  at  r  =  0: 

u(uu  t)  =  {  (uf  +  emA)  exp(— 3r )  -  emA}1/3  >  0. 

Note  that  the  proton  speed  is  reduced  significantly  below  its  initial  value  u  <  Ui  on  a  time 
scale  of  the  order  of  the  (p-e)  collision  time  ~  l/i^.  The  angular  spread  is  determined  by 
the  magnitude  of  a(u).  In  the  limit  ao  — ♦  oo  and  u3  >  u3  >  emA, 

,  N  3u3  1 

“  T.(i-«8MS)‘ 

Hence,  for  A  <  1  and  u  ~  1,  a(u)  >  1  and  the  angular  spread  of  the  proton  beam  remains 
small. 

On  the  other  hand,  for  A  >  1  and  u  ~  1,  a(u)  <  1  and  the  above  solution  is  no  longer 
valid.  L  and  M  investigated  this  limit  and  observed  growth  in  the  transverse  velocity  spread 
of  the  proton  beam  [1]. 


B.  Electron  Drag  Balanced  by  Inductive  Electric  Field  E  ^  0  and  Ve  ^  0 


In  the  case  where  an  inductive  field  E  ^  0  (or  an  equivalent  “external”  force)  is  present 
to  balance  the  drag  forces  on  the  protons,  and  Ve  ^  0,  we  have 


dFp 

dr 


u2  du  \ 


q{u)  + 


F  .  p  ,  A  d  dFp 

2  )Fp  +  ud  dee  Fp  +  2v?e  dee  do 


+SoQ{t)5(u  -  Ui)2aoexi>(-ao02/2) 
q{u)  =  u3  —  v? 


(17) 

(18) 


where  r  =  upet,  u  —  v/vq,  and  vq  is  now  defined  to  be  the  proton  beam  speed  at  which  the 
electron  drag  is  balanced  by  the  inductive  electric  field  eE  =  mpUpe^o  —  Ve).  We  assume 
A  <  1  and  we  ignore  the  boron  drag  emA  ;$  1. 
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For  values  of  u  such  that  q(u)  >  (u9)2/2,  the  solution  for  Fp(u,  r)  is: 


Fp(u,t )  =  Fp(u,  r)2a(u)  exp{— a{u)92/2) 


(19) 


where 


Fp{u,  r)  =  {G(u  -  u{uu  t))  -  G(u  -  Ui)} 


q(u) 


a(u)  = 


a0h2(u ) 


(1  +  a,oAhi(u)) 


hi(u)=  f 

Ju 


uq{u) 

h2  (u)  =  exp  (2  [  du-y-r  ) 

V  QWj 

and  u{ui,r)  is  the  solution  of  the  characteristic  equation  ^  ^  with  initial  conditions 

u  =  Ui  at  r  =  0: 

u(v,i,  r)  =  1  +  {ui  —  1)  exp(-r). 

The  proton  velocity  u  decreases  from  it*  and  approaches  u  — »  1. 

In  the  limit  do  — *•  oo  and  1  >  (u  —  1)  >  92/ 2,  we  have  a(u)  ~  2/A  >  1,  and  the  angular 
spread  remains  small. 

For  (u  —  1)  ~  02/ 2,  the  above  solution  is  no  longer  valid.  The  Fokker-Planck  equation 
for  a  proton  distribution  sharply  peaked  about  u  =  1  can  then  be  approximated  by 


SFp  _  S_  (  f)  1  8,2  1  £,8^. 

8t  85r+2/  p+eaee  p+eaee9»‘ 


(20) 


where  we  have  introduced  new  dimensionless  variables,  6  =  (2/ A) 1^29  and  £  defined  by 
u-l  =  (A/2)(. 

We  do  not  have  an  analytic  solution  of  this  equation.  However,  the  averaged  distribution 
function  Fp  =  /  d£Fp,  is  determined  by: 


OF, 

dr 


=  lS=tF 

Odd6  p+6de6d9  J 


and  the  time  independent  solution  of  Fp  is: 


FP  — ►  exp(—92/2)  =  exp(—92/A). 
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The  transverse  velocity  spread  v\  of  the  proton  beam  is  therefore: 

~  \ 

A  v. 


pe 


A  similar  result  is  also  obtained  from  a  consideration  of  the  beam-envelope  equation  de¬ 
rived  in  the  next  section  with  the  inclusion  of  a  field-reversed  magnetic  field  and  a  transverse 
electrostatic  field. 


C.  Beam-envelope  Equation 


By  taking  appropriate  moments  of  the  Fokker-Planck  equation,  we  derive  a  set  of  coupled 
equations  to  describe  the  mean  square  width  of  the  proton  beam  [5]. 

Let  the  mean  square  width  of  the  proton  beam  X 2  be 

f  d3r  J  d3vFx2 


X2  =  (x2)  = 


J  d3r  /  d3vF  ' 

Then,  from  the  Fokker-Planck  equation  (see  Appendix  B),  we  have: 


(21) 


Jt(  t)-<ot*>=° 


d_ 

dt 


(xvx)  -  (vl)  -  (xEx  -  — =  {xvxcpe)  +  { XVxCpB ) 

It  (”■) "  (l  (VlEl  -  ^T~))  -  (v?c’*)  + 

(tf)  =  (Vl°pe)  +  ( VlCvB ) 


C 

dt 


(22) 

(23) 

(24) 

(25) 


We  now  make  the  ansatz  that  the  x-component  of  the  particle  velocity  vx  be  expressed 
in  the  form 

x  dX 


Vx  = 


+  6vx 


x  X  dt 

and  we  assume  that  Ex  and  Bv  vary  linearly  with  x: 


(26) 


Then  from  Eq  (22), 


Ex  =  E'xx  By  =  B'x 


(xSvx)  =  0 
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and  we  have 


{xv.)=Xd-£ 

<*>-(£V+M 


(£ 

\m 


( -  (xEx 
(vxEx  - 

where  the  betatron  frequency  ftp  is 


\dt  J 

xv  ZB. 


*)>- 


2  V2 


QiX 


(- 

\m 


S¥l))—BixW 


nl=^-yp-—K 


mpc 


and  the  mean  velocity  of  the  proton  beam  is  Vp  =  Vpz. 

Prom  Eq  (23)  and  Eq  (24),  we  have 

J2  v 

X-p-  +  fi^X2  =  (( Svx )2)  +  (xvxCpe)  +  ( xvxCpB ) 


=  ((SVX)2)  -  (fpe  +  VpB)X—  + 


dt 


dX 


(27) 


(28) 


.dX 


|*2  ((fo,)2)  =  X2  (vlc„)  +  X2  (vICpB)  -  2  (iVpCp.)  X—  -  2  (xv,CtB)  X— 


dt 


=  -2 upeX2  {(Svx)2)  +  vpBX‘ 


mB 


V2  + 

mp  +  mB  v 


(29) 


Jt  ((%)2)  =  —2Upe  (K)2)  +  ^X5 


™fi . Vp2  + 

mp  +  mB  p 


(30) 


where  we  have  assumed  that  Vp  >  {((Svx)2) ,  ((uy)2)}  >  Te/mp. 

The  square  root  of  the  mean  square  equilibrium  width  X,  neglecting  collisions,  is: 

<(fe.)2)1/2 


X  = 


Sip 


In  the  presence  of  collisions,  the  proton  beam  expands  due  to  p-B  collisions.  However,  a 
stationary  state  can  be  achieved  if  the  heating  of  the  proton  beam  due  to  p-B  collisions  is 
balanced  by  cooling  due  to  electron  drag.  At  steady  state,  we  have: 


((<M2)  +  (ivy)2)  _  vl  vvB  _  , 

. — .1  ■■■■  11  1  * _ _  ~  _  I 

w  rO  T  * 


V 2 

p 


V 2 
v 


3(10)* 


pe 


/  rp  \  3/2 

|  —L  I  Ub 
\£pB  / 
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In  summary,  the  transverse  velocity  spread  ^  of  the  proton  beam  remains  small  if 
_2 

i/ps/^pe  <  1.  For  ~  VpB/vpe  <  0.25,  the  electron  temperature  and  the  density  ratio  of 
borons  to  electrons  are  bounded  by  the  inequality 

(^)3/2—  <1.9(10)-4. 

\t-pB  /  Tle 

IV.  SUMMARY  AND  CONCLUSIONS 

In  the  proton-boron  colliding  beam  fusion  reactor  design  concept  of  Rostoker,  Binder- 
bauer,  and  Monkhorst,  proton  and  boron  beams  collide  with  a  relative  resonant  energy  Spb 
of  collision  equal  to  Zvb  ~  600  KeV.  The  fusion  reactivity  decreases  if  the  temperatures  of 
the  fuel  ion  beams  acquire  a  thermal  spread  or  if  the  interaction  energy  deviates  from  the  res¬ 
onant  energy.  Efficient  operation  therefore  requires  that  optimal  conditions  be  maintained 
for  times  long  enough  for  a  significant  number  of  fusion  reactions  to  occur. 

Quasi-stationary  colliding  beam  configurations  can  in  principle  be  sustained  for  times 
longer  than  Coulomb  collisional  times.  Collisional  drag  forces  reduce  the  relative  velocity 
between  particle  species.  However,  if  the  drag  forces  are  balanced  with  inductive  electric 
fields  and/or  “external”  forces  acting  on  the  proton  and  boron  beams,  the  energy  of  collision 
of  the  fuel  ion  beams  can  be  maintained  close  to  the  resonant  energy.  It  is  also  essential  to 
control  the  transverse  velocity  spread  of  the  proton  beam  induced  by  proton-boron  collisions. 
This  can  be  achieved  by  means  of  electron  drag  if  the  boron  density  is  sufficiently  low. 

However,  in  order  to  sustain  steady  state,  energy  must  be  supplied  at  a  rate  sufficient 
to  replace  the  streaming  energy  transformed  into  thermal  energy  by  particle  collisions.  This 
energy  supply  constitutes  circulating  power  in  a  reactor. 

From  the  momentum  balance  equations,  the  magnitude  of  the  inductive  electric  fields  or 
the  “external”  forces  necessary  for  steady  state  can  be  calculated.  The  power  which  must 
be  supplied  can  then  be  readily  estimated. 

We  find  that  in  the  case  of  inductive  electric  fields,  the  ratio  of  the  circulating  power 
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Pc  =  Pd  (required  to  maintain  the  quasi-stationary  state)  to  the  fusion  power  Pf  is 

Pc  f  ne  ( Z  +  Sl  7ieZ(l  +  S)2  1 

P/  ‘  \nBZS(Z-l)2  +  npS(Z-l)2  ^  j 

where  the  parameter  S  is  defined  to  be  S  =  5^7*7  ~  264  Te  is  the  electron  tem¬ 

perature,  Z  the  boron  atomic  number,  n,  (i  =  e,p,  B)  the  species  density,  m,  (i  —  p,  B)  the 
species  mass,  vpb  the  proton-boron  collision  frequency,  and  Upe  the  proton-electron  colhsion 
frequency. 

The  absolute  minimum  value  of  Pc/P/  is  35,  which  corresponds  to  a  density  ratio  of 
nB/ne  —  0.1  and  S  =  Z1//2  (or  equivalently  Te  =  25KeV).  The  large  magnitude  of  Pc/Pf  is 
undesirable  since  it  implies  the  need  for  very  efficient  means  of  converting  the  available  power 
into  electrical  power.  The  power  available  for  conversion  into  electrical  power  is  Pj  +  Pc. 
Then,  introducing  an  effective  conversion  efficiency  parameter  denoted  by  77,  we  find  that 
net  electrical  power  is  available  only  if  77 (Pf  +  Pc)  is  greater  than  the  power  Pc  required  to 
maintain  steady  state.  Thus  the  effective  conversion  efficiency  must  be  greater  than  97%, 
that  is,  77  >  35/36  ~  0.97. 

In  the  case  of  “external”  forces,  we  find  that  the  power  required  to  maintain  the  quasi- 
stationary  state  is 

This  power  is  lower  than  that  required  for  a  quasi-stationary  state  with  inductive  electric 
fields.  At  an  electron  temperature  of  Te  =  25KeV  and  for  the  optimal  density  ratio  of 
nB/ne  =  0.1  (at  which  value  the  fusion  power  is  maximized),  the  circulating  power  is 
a  factor  of  ~  9  times  the  fusion  power  {Pc/Pf  ~  9).  At  higher  electron  temperatures, 
the  magnitude  of  Pc/P/  is  reduced  and  tends  to  a  minimum  of  Pc/P/  — ►  5.1  in  the  limit. 
Te  >  42.6/(l  +  20nB/ne)2/3  KeV  (or  >  i )•  The  required  effective  efficiency  of 
conversion  to  electrical  power  must  therefore  be  greater  than  84%,  that  is,  77  >  5.1/6. 1  ~ 
0.84. 

Unfortunately,  at  the  high  electron  temperatures  necessary  to  keep  the  power  ratio  Pc/P/ 
close  to  its  minimum  value,  electron  drag  will  not  be  strong  enough  to  cool  the  proton 


tH1  1+ 


_ ZUe _ 

S(ne  +  Z(Z  —  1)715) 
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beam  (i.e.,  suppress  the  transverse  velocity  spread)  unless  the  boron  density  is  considerably 
below  its  optimal  value  of  715 /ne  =  0.1.  From  steady-state  solutions  of  the  Fokker-Planck 
equation  for  the  proton  distribution  function,  we  estimate  that  the  transverse  proton  beam 


temperature  T±  is  given  by 


Zn*s  _  1-32(10)»2« 

&pB  Vpe  Tle  Tig 


At  Te  ~  68KeV,  electron  drag  will  be  effective  in  limiting  the  transverse  proton  beam 
temperature  {T±/£pb  <  0.25)  only  if  the  density  ratio  is  <  0.005.  At  this  low  boron 

density  of  riB/ne  =  0.005,  the  fusion  power  is  reduced  to  10%  of  the  maximum  possible.  The 
power  ratio  is  Pc/Pf  ~  7.4  at  Te  =  68KeV  and  tib/t^  =  0.005,  and  the  required  effective 
conversion  efficiency  would  then  have  to  be  greater  than  88%,  that  is,  77  >  7.4/8.4  =  0.88. 

The  preferred  scenario  is  the  one  involving  “external”  forces.  It  will  be  essential  to  have 
a  highly  efficient  means  of  generating  the  “external”  forces  necessary  to  transfer  momentum 
to  the  ion  beams.  It  is  as  yet  unclear  how  this  can  possibly  be  accomplished.  But  even 
if  the  means  for  efficiently  generating  the  “external”  forces  were  available,  the  ratio  of  the 
circulating  power  to  the  fusion  power  is  still  very  high,  Pc/Pf  >5.1,  which  implies  that 
effective  conversion  efficiencies  greater  than  84%  are  required. 

It  should  be  noted  that  even  this  high  value  of  77  ~  0.84  still  underestimates  the  required 
effective  conversion  efficiency  for  the  following  reasons: 

(1)  In  our  simplified  calculation  of  77,  we  assumed  that  the  mechanism  for  transferring 
momentum  to  the  proton  and  boron  beams  to  make  up  the  losses  due  to  collisions!  drag  is 
100%  efficient.  It  is  unlikely  that  this  level  of  efficiency  can  ever  be  achieved. 

(2)  We  also  assumed  that  the  conversion  efficiency  is  the  same  uniformly  high  value  for 
all  the  various  forms  of  the  energy  available  for  conversion  into  electricity  (e.g.  thermal 
energy,  directed  energy,  bremsstralilung).  If  the  transverse  spread  of  the  proton  beam  is 
limited  by  electron  cooling,  the  dissipated  energy  is  eventually  transferred  to  electron  energy 
and  must  be  removed  through  the  electron  loss  channels.  At  the  above  nominal  electron 
temperatures  and  boron  densities,  the  bremsstrahlung  loss  rate  is  much  smaller  than  the 
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dissipated  power,  and  a  large  fraction  of  the  electron  energy  will  have  to  be  removed  as 
electron  thermal  energy.  The  conversion  efficiency  for  thermal  energy  is  ~  0.4,  well  below 
the  required  effective  conversion  efficiency  of  rj  >  0.84.  At  the  present  time,  direct  energy 
converters  (which  convert  energetic  charged  particle  into  electricity)  have  been  designed  with 
conversion  efficiencies  no  higher  than  ~  75%. 

In  addition,  PJPf  ~  5.1  is  possible  only  at  high  electron  temperatures.  However,  at  such 
temperatures,  the  boron  density  will  have  to  be  well  below  its  optimal  value  ( TiB/ne  <  0.1) 
if  electron  drag  is  to  be  effective  in  limiting  the  transverse  proton  beam  temperature.  The 
fusion  power  Pj  is  then  only  a  fraction  of  the  maximum  possible. 

In  summary,  we  conclude  that  the  proton-boron  colliding  beam  fusion  reactor  is  not  a 
viable  concept  unless  technology  capable  of  very  high  energy  conversion  efficiencies  (no  less 
than  84%)  can  be  developed. 
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Appendix  A:  Mixture  of  Inductive  Electric  Field  and  Non-inductive  “External” 

Forces 


We  now  consider  a  “steady  state”  maintained  by  the  presence  of  an  inductive  electric 
field  E  of  arbitrary  magnitude  and  by  additional  “external”  forces  forces  Tp  and  TB  which 
act  to  maintain  the  streaming  velocities  of  the  ions  at  fixed  prescribed  values.  We  take  E 
to  be  a  parameter  and  we  establish  that  the  dissipated  power  is  minimized  when  E  =  0. 

The  modified  momentum  equations  are: 


71.B 


Tp  +  eE  =  rripUpeCVp  -  Ve)  +  mBuBp(Vp  -  VB) — 


n„ 


Tb  +  ZeE  =  mpZ2Upe(VB  —  Ve )  +  mBuBp(VB  —  Vp ) 


-eE  =  mvVpt{Ve  -  Vp)^  +  mpZ\e(Ve  -  VB )— 


nB 


n{ 


(Al) 

(A2) 

(A3) 


Solving  for  (Vp  -  Ve),  (VB  -  Ve),  TV,TB  in  terms  of  Vp  -  VB  and  E,  we  obtain 
vp  -  K  =  7  ,  (z2nB(Vp  -  VB)  +  ^5- 

( np  +  Z2nB)\  H  rnpvpej 


VB-Ve  = 


(' np  +  Z2nB)  \ 
Z2nB 


np(Vp  -  VB)  ~ 


eEne ' 
mpupe/ 


{  np 


.  +  PnB 


Vb)~ 


eEuBZjZ  -  1) 

np  +  Z2nB 


= 


Tin 


(A4) 

(A5) 

(A6) 


The  dissipated  power  is 


Pd  =  ripFpVp  +  nBFBVB  +  eE{npVp  +  ZnBVB  -  neVe) 


=  mBnBuBp(Vp  -  VB )2 


1  + 


Zne 


S(np  +  Z2nB) 


+ 


(( sEne)2 


TUpUpefop  +  Z2nB) 


(A7) 


Note  that  finite  values  of  E  increases  the  magnitude  of  PB- 
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Appendix  B:  Fokker-Planck  Equation 


The  Fokker-Planck  which  determines  the  time  evolution  of  the  distribution  function 
Fa(r,v,t)  of  the  particle  species  in  a  collisional  plasma  is  [3]: 


+  - -(e+ 

at  or  ma  \ 


vxB\  w.j.cwa.a, 


(Bl) 


where  the  collision  operator  Cab  for  multiple  small  angle  Coulomb  scattering  between  species 
“a”  and  species  “b”  is  given  by: 


and 


Cab  —  Fab 


(i  .  IM  A  (F  + 1  dl 

V  mb)  dva  \  adva)  2dvadvp 


r ab  — 

Ane2aelXab 

o 

ml 

hb  — 

f  A/ 

J  |v  —  v'\ 

9b  = 

f  d3v'  Fb(v')\v 

The  subscripts  a,  /?  of  va,  vp  denote  the  cartesian  components  x,  y,  z  of  the  velocity  vector 
v.  A  sum  over  repeated  indices  a,  f3  is  implied.  ea  and  ma  are  the  charge  and  mass  of  species 
“a”  and  Xab  the  Coulomb  logarithm. 

Particle  number,  total  momenta,  and  total  energy  are  conserved  during  collisions, 
f  d3vCab  =  0,  /  d3vmavCab  +  Jd3vmbvCba  =  0,  \  J  d3vmav2Cab  +  \  /  d3vmbV2Cba  —  0. 

For  proton-electron  collisions,  where  Fe  is  taken  to  be  a  drifting  Maxwellian  Fe(v)  = 
ne  (^r)3/2exp(— me(v  —  Ve)2/2 Te)  and  the  electron  thermal  velocity  is  much  larger  than 
typical  proton  velocities  i t(me/2Te)^2  ;$  1,  we  can  approximate  the  collision  operator  Cpe 
by: 


r  -V 

^pe  —  1  pe  g 


©"©)"{» •S's'-.- 

f  d  „  Te  d2  „  1 
Upe  {  dva  pUa  +  mp  dvQdva  p) 


+ 


me  dv. 


5L_f\ 

adva  Pj 


+ 
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where  ua  =  va  -  V*  and  vpe  =  Tpe^  (|)1/2  (^)3/2  ^  is  the  proton-electron  frequency. 

For  proton-boron  collisions,  where  Fb  is  a  drifting  Maxwellian  and  the  boron  thermal 
velocity  is  much  less  than  typical  proton  velocities  vims /IT b)1^  1,  we  have: 


C„r  =  TvbNe 


=  VpB  ^0 


mB  +  mp 


(i + I + ■  ■  • 

ms  dva  v?  2dvadv$  u  \  u2  ) ) 

fmpfl  ttq  lfll/.  UgUp\  dFp) 

\mBdva  pu3  2  dva  u  V  a<3  v?  )  dvp  ] 


where  ua  =  vQ-  VB,  vpB  =  TpBNB{l  +  mp/mB)/vl  is  the  proton-boron  collision  frequency, 
and  vq  is  the  characteristic  speed  \VP  —  VB\  of  the  protons  relative  to  the  borons. 
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FIGURE  CAPTIONS 


FIG.  1.  Quasi-stationary  state  maintained  by  inductive  electric  field;  Pd/P{  vs  Te(KeV) 
for  rib/ne  =  0.05,0.1,0.15. 

FIG.  2.  Quasi-stationary  state  maintained  by  non-inductive  “external”  forces;  Pd/Pj  vs 
Te(KeV)  for  nb/ne  =0.01,0.1,0.15. 
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(1)  nB/ne=0.05 

(2)  nB/ne=0.10 

(3)  nB/ne=0.15 


Figure  2 
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Abstract:  MHD  stability  of  highly  elongated  tokamaks  (termed  a  belt  pinch)  are  considered  for  high 
bootstrap  fraction  cases.  By  employing  high  triangularity  or  indentation,  and  invoking  wall 
stabilization,  and  P  can  be  increased  by  a  factor  of  roughly  3  by  increasing  k  from  2  to  4. 
Axisymmetric  stability  up  to  K  =  4  tolerable  by  employing  a  shell  which  conforms  more  closely  to  the 
boundary  than  in  present  experiments.  Engineering  difficulties  with  a  close  fitting  shell  in  a  reactor 
environment  may  be  overcome  by  employing  a  liquid  lithium  alloy  shell.  Rapid  metal  flows  can  lead 
to  potentially  deleterious  plasma  shifts  and  damping  of  the  flow. 


Several  authors  have  found  that  kink  stability,  as  measured  by  PN,  decreases 
for  k  >  2  [1,2].  We  have  used  the  MHD  stability  code  PEST  [3]  together  with  the 
equilibrium  code  JSOLVER  [4]  to  investigate  kink  stability.  Unlike  previous 
investigations,  we  examine  high  bootstrap  fraction  profiles  with  greater  than  99% 
pressure  driven  current..  The  plasma  shape  is  parameterized  by  k  and  8  [4].  We  have 
also  considered  indented  cases;  the  boundary  shape  is  specified  bending  a  purely 
elliptical  plasma  inward  along  a  constant  radius  to  produce  the  desired  indentation  A.. 
A  wall  is  assumed  at  b/a  =  1.3.  The  results  can  roughly  be  summarized  by  stating  that 
with  sufficient  triangularity  cases  with  k  =  3  have  stable  (3  values  twice  as  high  as 
similar  cases  for  k  =  2.  If  there  is  sufficient  triangularity,  cases  with  k  =  4  have  three 
time  higher  stable  (3.  Cases  with  k  =  4  are  very  challenging  numerically  and  it  has 
only  been  possible  to  examine  stability  up  to  mode  numbers  n  =  5. 


In  the  following,  we  have  required  that  modes  with  n  =  5  be  stable  without  a 
wall.  This  is  a  plausible  restriction,  since  the  use  of  wall  stabilization  implies  that 
resistive  wall  kink  modes  must  be  stabilized.  In  a  reactor  environment,  plasma 
rotation  is  not  a  viable  stabilization  method  because  the  required  rotation  velocity  is 
too  high.  Feedback  stabilization  is  under  examination  for  low  n  modes.  The  unstable 
eigenfunctions  for  the  high  n  cases  extend  far  into  the  plasma,  indicating  that  they  are 
not  edge  localized  modes,  but  rather  would  have  serious  consequences.  Thus  they 
would  require  feedback  stabilization.  Practical  considerations  probably  will  restrict 
the  number  of  n  modes  which  can  be  stabilized.  Thus,  we  have  determined  the  stable 
beta  values  for  cases  where  modes  with  n  =3  and  above  are  stable  without  a  wall,  and 
thus  would  not  require  feedback.  For  triangularity  =  .71,  resulting  stable  (3*  for  k  = 


2,3  and  4  is  4.4%,  7.2%  and  8%,  respectively,  and  is  10%  for  a  case  with  indentation 
=  .5  and  k  =  4.  Requiring  (with  no  wall)  for  modes  with  n  =5,  the  stable  P*  for  k  = 
2,3  and  4  is  6.6%,  10.5%  and  14%,  respectively,  and  is  also  14%  for  the  indented  case 
with  k  =  4.  Raising  triangularity  even  slightly  greatly  improves  the  high  elongation 
cases;  for  indentation  .78,  the  stable  beta  values  for  for  k  =  2  and  3  are  6.7  %  and 
13%,  respectively.  For  k  =  4,  the  stable  P*  increases  to  22%. 

At  high  elongation,  axi-symmetric  instabilities  require  stabilization  by  a 
passive  shell,  and  an  active  feedback  system.  For  active  stabilization  to  be  feasible 
the  resistive  growth  rate  cannot  be  too  large.  An  initial  value  code  has  been  written  to 
examine  resistive  wall  axi-symmetric  instabilities.  Since  the  resistive  wall  evolution 
is  slow  compared  to  the  Alfven  time,  the  perturbed  Grad-Shafranov  equation  is  used; 
for  the  stability  calculation,  we  linearize  in  the  perturbed  flux,  starting  from  a  full 
numerical  equilibrium  from  the  code  TOQ.  The  plasma,  wall  and  external  coils  are 
broken  into  1000  -  2000  finite  elements,  so  that  all  plasma  deformations,  wall 
responses  and  feedback  effects  can  be  described  accurately.  Convergence  in  the 
number  of  elements  was  checked  routinely  in  the  following  calculations. 


n=0  Resistive  Wall  Growth  Rate  vs.  Elongation 


Fig.  1:  Normalized  growth  rate  vs.  plasma  elongation  ( triangularity  .7 )  for 
conformal  and  rectangular  shells.  The  minimum  distance  from  the  plasma  is.l  a 
on  each  side  of  the  rectangular  shell,  and  as  indicated  for  th  e  conformal  shells. 
The  plasma  l,  -.7,  and  fip  =  0,  or  fip  =  2.. 5  if  indicated. 


Growth  rates  normalized  to  the  shell  L/R  decay  time  (for  anti-symmetric  decay)  are 
shown  in  figure  1,  for  conformal  shells  and  a  rectangular  shell,  for  a  plasma  with  8  = 
.7.  The  tokamak  facilities  DIH-D  and  TDV,  which  have  explored  high  elongation, 
have  vacuum  vessels  which  can  be  regarded  as  roughly  rectangular,  for  highly 
triangular  plasmas.  As  can  be  seen,  a  close  conforming  shell  enables  much  higher 


elongation  without  high  growth  rates,  as  compared  to  a  “close”  rectangular  shell. 
Also  note  that  advanced  tokamak  operation  (  Pp  ~  2.5  )  further  enables  high 
elongation. 

For  a  conducting  shell  corresponding  to  2  cm  of  lithium  with  reactor 
dimensions,  Twall  -  .1  sec,  so  the  growth  rates  low  enough  to  be  stabilized.  We  have 
found  positions  for  sensor  and  active  loops  which  give  stabilization,  for  two  or  more 
sensor  and  active  coils  located  at  the  top  of  the  vessel,  (and  two  at  the  bottom).  The 
active  coils  are  located  1.3  meters  away  from  the  plasma,  far  enough  to  be  behind  the 
shield.  Stability  is  sensitive  to  the  size  of  the  hole  in  the  shell  which  is  needed  for  a 
divertor;  hole  widths  ~  .7  are  acceptable  for  d  =  .7,  but  not  much  beyond  this. 
Indented  cases  have  much  lower  wall  growth  rates,  and  less  sensitivity  to  the  divertor 
hole. 


In  a  reactor  environment,  there  are  serious  engineering  difficulties  in 
designing  a  solid  conducting  shell  close  to  the  first  wall.  Even  a  few  cm  of  material 
can  seriously  degrade  the  tritium  breeding  ratio,  and  mechanical  stress  from  neutron 
heating  and  embrittlement  can  lead  to  low  reliability.  In  reactor  design  studies  such 
as  ARIES  RS  [8],  the  conducting  shell  is  placed  partially  or  completely  behind  the 
blanket,  limiting  the  elongation  and  beta.  Liquid  lithium  alloys  used  as  a  conductor 
can  conceptually  solve  these  problems.  Several  liquids  have  breeding  ratios  greater 
than  one,  and  are  not  subject  to  mechanical  stress  concerns. 


Here  we  will  limit  ourselves  to  the  question  of  whether  a  liquid  conductor  in  a 
magnetic  field  can  act  as  a  stabilizing  shell.  We  consider  a  thin  shell  of  fluid  with 
finite  conductivity,  soE  +  vxB  =  pj.  The  fluid  could  either  have  solid  walls  on  both 
sides  (i.e.,  be  in  channels  just  behind  the  first  wall)  or  have  a  free  surface  facing  the 
plasma.  We  use  low  aspect  ratio  reduced  MHD  to  describe  the  liquid,  and  consider 
the  coupling  of  the  liquid  to  plasma  MHD  perturbations.  In  response  to  an  external 
time  varying  magnetic  field  of  frequency  to,  the  liquid  attempts  to  respond  by  flowing 
with  the  inductive  E  x  B  velocity.  Note  that  in  such  a  case  the  eddy  currents  (which 
tend  to  stabilize  the  plasma)  which  would  arise  in  a  solid  wall  would  not  be  produced, 
since  E  would  be  balanced  by  v  x  B  in  the  case  of  a  liquid,  rather  than  of  p  j  as  in  a 
solid.  However,  the  liquid  cannot  flow  into  the  wall,  which  results  in  a  boundary 
layer  of  width  (topp/B2)  /  km  (where  k,u  is  the  parallel  wave  number,  and  all  quantities 
are  in  MKS  units).  Within  this  layer,  v  x  B  cannot  cancel  E,  so  eddy  currents  are 
produced  as  in  a  solid.  A  shell  of  fluid  which  is  thinner  than  this  behaves  as  a  solid. 
In  the  case  of  the  vertical  instability,  co  is  several  tens  of  Hz,  and  k|,for  the  lowest 
Fourier  Harmonic  is  ~  q  R.  For  lithium,  this  width  is  -  10  cm.  For  Li  alloys  (which 
are  denser)  this  width  is  thicker.  Thus,  even  for  the  higher  harmonics  in  the  system,  a 
2  cm  liquid  layer  should  behave  as  a  solid  shell.  We  expect  that  other  mechanical 
obstructions  to  the  flow  (such  as  channels)  reinforce  this  conclusion.  Space  does  not 
permit  a  discussion  of  other  effects  such  as  rapid  liquid  flow,  and  ways  to  reduce  the 
flow  damping  and  plasma  shifting  to  low  levels  in  fast  flows  to  acceptable  levels. 
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II.  Effects  of  Liquid  Metal  Walls  on  Equilibrium  and  Stability 

A.Y.  AYDEMIR,  Institute  for  Fusion  Studies 

For  a  flowing  liquid  metal  wall  to  be  effective  in  controlling  MHD 
instabilities,  the  flow  needs  to  satisfy  two  essential  requirements.  The  first  is  an 
equilibrium  condition:  In  steady  state,  the  flow  has  to  be  of  the  form  u  =  F  (\|t)  B/p  + 
R2£2(i|/)V£  where  F  (vy)  and  £2(\j/)  are  arbitrary  flux  functions1,  and  £  is  the  toroidal 
angle.  The  flow  is  within  flux  surfaces.  It  is  also  easily  seen  that  u  can  be  purely 
toroidal  but  cannot  be  purely  poloidal. 

The  second  is  a  stability  requirement:  For  the  liquid  metal  wall  to  respond  as  a 
"perfect"  conductor  to  MHD  instabilities,  the  flow  has  to  satisfy  QR/  £2W>  1,  where 
£>Ris  an  effective  rotation  frequency  for  the  liquid  wall,  and  £2W  =  2n  /  Tw,  and  Tw  is 
the  resistive  wall  time  on  which  the  resistive  wall  modes  are  expected  to  grow. 
Satisfying  the  equilibrium  and  stability  conditions  simultaneously  is  a  nontrivial  task 
that  does  not  seem  to  have  received  adequate  attention.  Violating  the  equilibrium 
condition  while  trying  to  satisfy  the  stability  requirement  will  drive  very  large 
currents  in  the  liquid  metal,  since  the  flow  will  drag  along  any  "frozen-in"  ambient 
equilibrium  field  (poloidal  or  toroidal).  Consequences  of  violating  the  equilibrium 
condition  are  graphically  demonstrated  in  the  figure  below.  The  pressure  contours  in 
an  elongated,  n=0  unstable  equilibrium,  depict  the  induced  rotation  in  the  plasma 
when  the  liquid  starts  rotating  poloidally  (but  not  within  flux  surfaces)  with  QR/ Qw= 
0.2,  due  to  induced  co-rotating  toroidal  currents  in  the  liquid.  Even  when  a  poloidal 
flow  is  almost  entirely  within  flux  surfaces,  it  leads  to  large  poloidal  currents  in  the 
liquid  wall,  Jp  >  10  Jeq  for  QR/ Qw=  1  due  to  dragging  of  the  toroidal  field.  These 
currents  would  be  expected  to  interfere  and  stop  the  flow. 
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